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We introduce a general (k,l) clique community, which consists of adjacent k-cliques sharing at least l ver- 
tices with k — 1 > l > 1. The emergence of a giant (k,l) clique community indicates a (k,l) clique 
percolation, which is studied by the largest size gap A of the largest clique community during network 


evolution and the corresponding evolution step Te. 


For a clique percolation, the averages of A and Te 


and the root-mean-squares of their fluctuations have power law finite-size effects whose exponents are re- 
lated to the critical exponents. The fluctuation distribution functions of A and Te follow a finite-size scal- 
ing form. In the evolution of the Erdés-Rényi network, there are a series of (k,l) clique percolation with 


(k,l) = (2 


, 1), (3, 1), (3, 2), (4, 1), (4, 2), (5, 1), (4, 3), and so on. The critical exponents of clique percolation 


depend on J, but are independent of k. The universality class of a (k, l) clique percolation is characterized alone 


by l. 
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Community structure is of great interest in the studies of 
networks . The term network community is defined as a 
group of vertices that are more densely connected each other 
than other vertices in a network. The clique and core fal 
are the two examples of network community. To analyze the 
overlapping community structure of networks, Palla et al. Ís) 
proposed the clique percolation method (CMP) to build up the 
communities from k-cliques, which is a fully connected sub- 
set of k vertices. Two k-cliques are considered to be adjacent 
if they share k — 1 vertices. A clique community is defined as 
the maximal union of k-cliques that can be reached from each 
through a series of adjacent k-cliques. 


When the size of a clique community is comparable to net- 
work size N, a giant clique community emerges and there is 
a clique percolation. Derényi et al. iel studied the clique per- 
colation of the Erdős-Rényi (ER) model|7]. The k-cliques 
percolation takes place when the probability of connecting 
two vertices in the network reaches the threshold pe(k) = 
[(k — 1)N]-MER-D [6]. The normal percolation transition 
of ER model corresponds to the k = 2 clique percolation. 


In this Letter, we introduce a general (k,l) clique commu- 
nity, where any two adjacent k-cliques share at least / vertices 
with k— 1 > l > 1. For the same network, the clique commu- 
nities of different (k, l) are different. As an illustration, (3, 1) 
and (3, 2) clique communities in a network are shown in Fig[I] 
and they are different. A (k, l) clique percolation appears with 
the emergence of a giant (k,l) clique community. The pre- 
vious clique percolation discussed in Ref.[6] corresponds to 
the special case 1 = k — 1 of general (k,l) clique percola- 
tion. The normal percolation transition corresponds to a (2, 1) 
clique percolation. Here we study the general (k,l) clique 
percolation of the ER model by analyzing the finite-size ef- 
fects of network evolution. From the power law exponents 
of finite-size effects, we can obtain the critical exponents of 
(k,l) clique percolation. 


At first, we calculate the threshold of the (k,l) clique per- 
colation using the approach developed by Newman, Watts, 


and Strogatz[9]. From the probability distribution pọ of ver- 
tex degrees, we define the generating functions Go(x) = 
Zro Pee” and Gi (x) = X po gea™, where qk = a 
1)/ < k > is the excess degree distribution of network 

The generating function Hı (x) of the probability distribu- 
tion for the sizes of the components reached from a randomly 
chosen edge satisfies a self-consistent equation 


Hı (x) = rG (Hı (2)). (1) 


When its first derivative H{ (1) becomes infinite, a giant com- 
ponent appears and there is a percolation. This corresponds to 
z =Gi (1) =1. 

After taking k-clique as the new unit and considering two 
cliques to be adjacent if they share at least / vertices, the aver- 
age degree of a k-cliques is 


z= 3 Cae ({;)2-©), (2) 


where p is the probability to connect two vertices with an 
edge. Correspondingly, the average excess degree of the k- 


FIG. 1: (color online) Sketches of different clique communities in a 
network. The vertices belonging to a clique community are marked 
with red color. The largest (3, 2) clique community is shown on the 
left and the largest (3, 1) clique community on the right. 
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FIG. 2: (color online) Log-log plot of the average size gap A and the 
root-mean-square of its fluctuation xa versus N for different (k, l) 
clique percolation. Power law behaviors of A and ya are confirmed 
by simulation data. From the slopes of fitting lines, 6; and 82 can be 
obtained and are summarized in Table I. 
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The threshold pe(k,l) of (k,l) clique percolation can be 
determined from the equation 


k-1 (2) 1 spri gk k" 
key aa [Ni | =1, 4) 


which gives 


k—l—1 


pelk, D = a(k, I) N= [1+ oF] , (5) 
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For 1 = k — 1, the correction term in EqJ5] vanishes and 
p(k, k — 1) = [(k - DN], which is in agreement with 
the result of Derényi et al id [1ah. For 1 < k — 1, the cor- 
rection term O(N ~ ET) exists. For general (k,l), only 
pelk, D) ~ N7 TF= is obtained by Bollobás and Riordan|[8}. 

In the following, we study (k,l) clique percolation of ER 
model with Monte Carlo simulation and finite-size effects of 
network evolution (14). During an evolution process, the 
largest (k, l) clique community has a size gap S1(T)—$1(T— 
1) at an evolution step T. The largest reduced size gap during 
the whole evolution process is 


1 
A = Smar ({S1(T) - S(T — 1)}), (7) 


which could be related to a percolation transition at the corre- 
sponding evolution step Te. From the results of A and Te in 


many Monte Carlo simulations, we can calculate the average 


size gap A and the average transition point Te. Since T = 
= 2 

p x N(N — 1)/2 and therefore Te x N°- TF=, it is conve- 

. . . = 2 

nient to introduce a reduced evolution step r = T/N?~ m= 


. . i, 2 
and the reduced transition point re = T/N 2-41, We an- 
ticipate that A and 7, have the power law finite-size effects 
as 


A(N) me, (8) 
Fe(N) = relo) ~ N74, (9) 


Using pe(k, l) in Eq. 5] we get 


rI (oo) = Ža(k,1). (10) 
The character of clique percolation is determined by the expo- 
nent (1. The clique percolation is continuous when 0 < 61 < 
1 and discontinuous when (3; = 0 (15). 

The fluctuations dA = A — A(N) and dr, = re — F-(N) 
are investigated also. Their root-mean-squares are defined as 


xa = Vv ([6A]?), (11) 
Xr = ([ôr.]2}, (12) 
which decay algebraically as 
a AN (13) 
Xr ~ N, (14) 


We anticipate a finite-size scaling form of fluctuation dis- 
tribution functions as 


Pa (ôA, N) 
P, (re, N) = 


NP fi(AN®?), (15) 
N2 falre N1”). (16) 


The universality class of continuous clique percolation is 
characterized by the critical exponents (1, 32, vı and v2. Dif- 
ferent clique percolation with the same critical exponents be- 
long to the same universality class. 

In Fig[2{a), the average A is plotted with respect to network 
size N for (k,l) = (3, 1), (3, 2), (4, 1), (4, 2), and (4,3). The 
log-log plot of A versus N show that A x N~®!, From 
the slope of fitting line, we can get the exponent 681. It has 
been obtained that 6; = 0.331(5) for (2,1), 6, = 0.32(2) for 
(3,1), 6, = 0.33(1) for (4,1), and 6, = 0.33(3) for (5, 1). 
Within error bars, 6; for l = 1 and different k are equal. At 
l = 2, we get 3; = 0.05(2) for (3,2) and 6, = 0.04(2) for 
(4,2). The exponents (3; of | = 2 are different from that of 
l = 1. Atl = 3, 6; = 0.01(1) for (4,3) and is different 
from that of | = 1 and l = 2. So the universality class of 
(k,l) clique percolation is characterized by l. In Fig[2{b), the 
log-log plot of xa versus N is shown for different (k,l). The 
slope of curve gives the exponent /2. It is found that {2 is 
equal to the corresponding (3, within error bars. 

The finite-size effects of 7.(N) — r-(co) are shown in 
Fig. Bla) for 1 < k — 1 and Fig. Bb) for] = k — 1. From 
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FIG. 3: (color online) Log-log plot of 7. — r-(oo) versus N. The 
results of | < k — 1 is shown in (a) and l = k — 1 in (b). The slopes 
of fitting lines give 1/1, which is defined in Eq.[9]and summarized 
in Table I. 
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FIG. 4: (color online) Log-log plot of xr versus N. The slopes of 
fitting lines give 1/v2, which is defined in Eq. [2]and summarized in 
Table I. 


the Monte Carlo data, the values of re(o0) and 1, are deter- 
mined simultaneously. For l = k — 1, we get 1/1, = 0.46(2) 
for (3,2) and 1/v, = 0.49(4) for (4,3). For] < k —1, 
we obtain 1/7; = 0.37(2) for (3,1), 1/1, = 0.38(2) for 
(4,1), and 1/1, = 0.67(20) for (4,2). The re(N) of (4, 2) 
clique percolation shows deviation from a simple power law 
of N. This deviation could be resulted by additional finite- 
size terms. 1/1, of different k and the same / are equal within 
error bars. We summarize r¢-(co) and 1/1 of different (k, L) 
clique percolation in Table I. 

The root-mean-square of dr, is shown in Fig/4] Our results 
of simulations confirm the power-law behavior of Xr. For k = 
3, we get 1/v2 = 0.33(5) for L = 1 and 1/v2 = 0.50(2) for 
l = 2. Atk = 4, we obtain 1/12 = 0.35(3) for l = 1, 1/v2 = 
0.52(4) for? = 2, and 1/v2 = 0.60(3) for l = 3. Other results 
of 1/2 are given in Table I. 1/12 of different k and the same 
l are equal within error bars. This confirms further that the 
universality of (k,l) clique percolation is characterized alone 
by l. The exponent vz agrees with vı at l = 1, 2, but differs 
from v at l = 3. 

For a network with size N, we simulate its evolution 
for many times. At each simulation, we can get the 
largest reduced size gap A and the corresponding transi- 
tion point re. From the results of many simulations, the 
fluctuation distribution functions Pa (dA, N) and P, (re, N) 
can be obtained. Five different network sizes N = 
100, 600, 1000, 6000, 10000 are chosen. Each fluctuation dis- 
tribution function has five different curves. Using the finite- 
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FIG. 5: (color online) Finite-size scaling functions of fluctuation 
distribution functions (a) Pa (A, N) and (b) P, (dre, N) for (4, 1) 
clique percolation. 


size scaling forms of Eq. [15]and[16] five curves of each dis- 
tribution function collapse into one curve of finite-size scaling 
function. The scaling variables are defined with the exponents 
Bg and v2 obtained from xa and xr. In Fig. [5] the finite-size 
scaling functions f;(SAN*?) and fo(ér.N!/”) are shown 
for (4, 1) clique percolation. fı and fz of other (k,l) clique 
percolation have similar behavior. 

In the evolution process of the ER model beginning from 
N isolated nodes to being fully connected with N(N — 
1)/2 vertices, there are a series of (k,l) clique percola- 
tion. The transition point of a (k,l) clique percolation is 
Te = a(k, DN EFi, which is of the same order of N 
for different clique percolation with equal k + 1. The se- 
quence of clique percolation in the ER model is as following: 
a) (2,1) clique percolation at Te = 0.5N, b) (3,1) clique 
percolation at T, = 0.5N4/3, c) (3,2) clique percolation at 
Te = 273/2 N3/2 ~ 0.354N3/2, d) (4,1) clique percolation 
at Te = 275/6 N3/2 ~ 0.561N3/?, e) (4,2) clique percola- 
tion at Te = 2712.571/5N8/5 ~ 0.416N8/5, f) (5, 1) clique 
percolation at T, = 27161/10N8/5 ~ 0.598N8/5, g) (4,3) 
clique percolation at Te = 271371/3N5/3 ~ 0.347 N5/3, and 
so on. 


For comparison with the investigations above, we investi- 
gate directly the size of the largest (k, l) clique community in 
a network to study general clique percolation. For a network 


TABLE I: 


Ta / N°7 TFT=T and critical exponents of (k,l) clique percolation. 
(oo) is obtained from Eq.D]and the analytic result r7 (o0) is given 
in Eq. E The results of e 1) clique percolation are taken from[14]. 


kD 1/v» 


(2,1) E 0.331(5) | 0.334(2) [0.331 (7) |0.334(3) 


Summary of reduced transition points re = 


(3,1)| 0.50(1) 0.32(2) | 0.36(4) | 0.37(2) | 0.33(5) 
(3,2)| 0.34(2) 2-3/2 | 0.05(2) | 0.03(2) | 0.46(2) | 0.502) 
(4,1)| 0.56(1) 275/6 | 0.33(1) | 0.35(2) | 0.38(2) | 0.35(3) 
(4,2)| 0.40(1) |2712.57'/5| 0.04(2) | 0.03(2) |0.67(20) | 0.52(4) 
(4,3)| 0.34(2) | 2-13-1/8 | 0.01(1) | 0.01(1) | 0.49(4) | 0.60) 
(5,1)| 0.60(1) | 2746*/2° | 0.33(3) | 0.38(4) | 0.30(5) | 0.37(6) 
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FIG. 6: (color online) Scaling size of the largest (4, 1) clique com- 


munity sı N?/” versus the reduced edge number r = T/N?’?. We 
take 6/v = 0.33. There is a fixed point at re = 0.559(4). 


with N vertices and T edges, the reduced size of the largest 
(k, L) clique community is denoted as sı (T, N). Near the tran- 
sition point Te of a (k,l) clique percolation, we anticipate a 
finite-size scaling form 
si(T, N) = NPS (tN), (17) 
where t = (T —Te)/Te and v is the critical exponent of corre- 
lation nes In the limit N > co, s1(T,oo) = OatT < Te 
iit E = A,t? at T > T. Our previous investigations 
ae mat the finite-size scaling form Eq. [7] for 
7 1) n percolation. 

At the transition point Te, s1ı(T, N)NP” |y Ln = 81(0) 
and is independent of network size N. Using T property, 
we can determine the transition point Te of (k,l) clique per- 
colation from the fixed point of sı (T, N)N°/”. 

In Fig. [6] we plot sı (T, N)N8/” of (4,1) clique commu- 
nity as a function of r for different N. 8/v = 0.33 has been 
taken. We find a fixed point at re = 0.559(4), which agrees 
with r.(00) = 0.56(1) obtained from network evolution. We 
expect that 8/v is equal to Bp (14). 

In summary, we introduce a (k,l) clique community con- 
sisting of adjacent k-cliques sharing at least | vertices with 
k—1 > l > 1. There is a (k,l) clique percolation with 


4 


the emergence of a giant (k,l) clique community. We study 
the (k,l) clique percolation by investigating the largest re- 
duced size gap A of the largest clique community during 
network evolution and the corresponding evolution step Te. 
If the average size gap A ~ N~®:, there is a continu- 
ous clique percolation for 0 < (6, < 1 and a discontin- 
uous clique percolation when 8; = 0. The reduced tran- 
sition point 7. obtained from Te has the finite-size effect 
Fe — reloo) ~ N~'/"1, The values of r.(0o) obtained from 
Monte Carlo data agree with the analytic result of Eq.[5] which 
is derived using generating function method. The sequence 
of (k,l) clique percolation in the ER model is as following: 
(2, 1), (3, 1), (8, 2), (4, 1), (4, 2), (5, 1), (4,3), and so on. The 
root-mean-squares of fluctuations 6A and dr, are found to de- 
cay algebraically as ya ~ N~* and yn ~ N7'/”. The 
fluctuation distribution functions follow the finite-size scal- 
ing forms Pa (ôA, N) = N®2 f,(SAN*2) and P,(dre, N) = 
N1/"2 fo(6r-N1/”), respectively. 

The universality of (k, l) clique percolation is characterized 
by the critical exponents (1, 82, 11 and v2. It has been found 
that the critical exponents of (k,l) clique percolation are in- 
dependent of k, but dependent on l. The universality class of 
(k,l) clique percolation is characterized alone by 1. 
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